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We study the role of the system-bath couphng for the generalized canonical thermalization [S. Popescu, et 
al., Nature Physics 2,754(2006) and S. Goldstein et al., Phys. Rev. Lett. 96, 050403(2006)] that reduces almost 
all the pure states of the "universe" [formed by a system 5* plus its surrounding heat bath B] to a canonical 
equilibrium state of S. We present an exactly solvable, but universal model for this kinematic thermalization 
with an explicit consideration about the energy shell deformation due to the interaction between S and B. 
By calculating the state numbers of the "universe" and its subsystems S and B in various deformed energy 
shells, it is found that, for the overwhelming majority of the "universe" states (they are entangled at least), the 
diagonal canonical typicality remains robust with respect to finite interactions between S and B. Particularly, 
the kinematic decoherence is utilized here to account for the vanishing of the off-diagonal elements of the 
reduced density matrix of S. It is pointed out that the non-vanishing off-diagonal elements due to the finiteness 
of bath and the stronger system-bath interaction might offer more novelties of the quantum thermalization. 
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Statistical mechanics is one of the most important and suc- 
cessful areas of modern physics. However its foundation is 
still debatable and is actually under debate. Most recently 
a mechanism for universal canonical thermalization has been 
found in the following meaning: almost all the pure states of 
the "universe" consisting of the considered system S and its 
surrounding heat bath B can be reduced into a generalized 
canonical state by tracing over the bath fl}]. Here, the allowed 
"universe" states lie in a subspace defined by a general con- 
straint R, which limits the number of the states to be finite, 
yet very large. It is noticed that in such generalized thermal- 
ization, the constraint R is rather general and need not be the 
energy shell of the "universe", and generally speaking it doses 
not lead to the usual canonical state in thermal equilibrium. 

When the constraint R is specialized as a total energy shell, 
the generalized canonical state becomes a usual canonical 
thermal state. Such a thermalization was described by the 
canonical typicality in ref. f?) associated with the so called 
overwhelming majority rule, based on the law of large num- 
bers, in counting the state numbers of the constrained "uni- 
verse" and its subsystems. Significant results concerning such 
canonical typicality have been obtained by several authors for 
different purposes Actually, to derive the 

canonical distribution of S from the microcanonical density 
matrix [9J or an entangled pure state [2J of the "universe" 
U = S + B, the interaction between 5' and B should be weak 
enough to allow a physical partition for the "universe". 

In this paper, we will quantitatively consider the effect 
of system-bath coupling on the above mentioned kinemat- 
ics of canonical thermalization leading to equilibrium canon- 
ical state. The present investigation only concerns the usual 
canonical state. In the weak interaction limit, the inverse tem- 
perature /3 — dS {E) /dE emerges from the thermodynamic 
entropy S {E) —hiVl {E, S) where SI {E, S) is the microstate 
number of the bath in the energy shell [E, E + S]. Since the 
interaction between S and B deforms the geometry of the en- 
ergy shell, we will reexamine the validity of the temperature 



definition from the new perspective offered by the general- 
ized thermalization |T, "2^] . We try to understand how the off- 
diagonal elements of the reduced density matrix of S vanish 
due to the factorization structure, in association with the ran- 
dom phase explanation jlOl ll ill . We find that if the mode 
number of the bath is not large enough and the system-bath 
coupling is strong enough, generally there exist nonzero off- 
diagonal elements in the reduced density matrix. Notice that 
these nonzero off-diagonal elements introduce quantum co- 
herence into the usual thermal equilibrium state and thus re- 
sult in novel thermodynamic features. 

We begin with a universal model: the system S we con- 
sider is an Af-level system with the Hamiltonian Hs = 
J2£n\n){n\, where \n) is the eigenstate with eigenvalue e„, 
n = 1,2,.., M; and the bath B is modeled as a collection of 
N harmonic oscillators of frequencies ojj {j — 1, 2, ..N) with 
the Hamiltonian Hb = J2j This model can be re- 

garded as a universal approach, because in the weak coupling 
limit, any heat bath could be universally modeled as a collec- 
tion of harmonic oscillators with the linear coupUngs to the 
surrounded system according to the proofs in ref. [13']. For 
this reason, the interaction Hj between B and S should be 
modeled to be linear with respect to the bath variables and 
aj. So we assume a simplest system-bath coupling 

Hi = Y,X^\n){n\{gja]+n.c.), (1) 

where A„ are real numbers. Note that Hi is of non-demolition 
since [iJg, Hj] = 0. Thus the interaction only causes the de- 
phasing of S, and the energy dissipation of S will not appear 

Obviously, the eigenvalues of the "universe" formed by B 
and S are 

JV 

£;(n, {nj}) = e„ (k) + (2) 
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FIG. 1: (Color Online) Geometrical deformation of the energy shell. 
Here we show explicitly the deformation, when S is a harmonic os- 
cillator (a) Without the interaction, the energy shell is the red area; 
(b) The interaction deforms the red area into the blue area. 

corresponding to the eigenstates |n, {nj}) = |n) (g) 

rijLi Wj ("■))' where \nj (n)} = D (ajn) \nj) is defined in 
terms of the Fock states \nj) of B and the coherent-state- 
generating operator D [ajn) = exp(aj„aj — h.c) with the 
displacement parameters ajn = ^KiQil (2ijJj); (k) = 
e„ — kA^j. Especially, the parameter k, — Y^j\gj^ / {Auj) 
reflects the role of the interaction between S and B. 

Note that the system-bath coupling deforms the energy 
shell of thickness 6 defined by the total constraint E < 
E {n, {uj}) < E+S. This energy shell determines a subset of 
the "universe" states. For convenience, we denote by V {E, 6) 
the subspace spanned by this subset. To have a clear picture of 
this deformation of energy shell, let us consider the following 
simple example: the system S* is a harmonic oscillator with 
frequency uj and A„ = n. In this case, the renormalized en- 
ergy of S is Es (n) = nuj — nn? . In Fig. [T] we illustrate the 
deformation of the energy shell when = 1. If there were no 
interaction, the above constraint would produce the red area. 
The interaction introduces the nonlinear term kv? to deform 
it into the blue area. 

Generally, in order to derive the canonical distribu- 
tion of S from an entangled pure state of the "uni- 
verse", we need to calculate the dimension JItv+i (S, 5, k) 
of V {E, S) and the dimension 0,^ [E — 6, k) of the sub- 
space {n, k) spanned by the states with the constraint 

fijUij < E + 5 — en (k) . (3) 
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Let us start with a pure entanglement state |?/'} = 
C (ri, {nj}) |ri, {nj}) of the constrained "universe", 
where denotes the summation under the total constraint. 
By tracing over the variables of B, the reduced density matrix 



Pg = Ttb dV') (V'I) of S' is obtained as 

Ps^Y^ Pn \n) {n\ + Fn,n \n) (m| , (4) 

where the diagonal elements are P„ = P„ (i?, n) — 
|C (n, nj)\ and stands for the summation under the 
constraint (|3]l. We represent this summation by the summation 
over the constrained indices [njjs „, then the off-diagonal el- 
ements can be written as 

Fn„.^ E E G("'".)^*K"^.)^™r4.)- 

where 

AT AT 

KSi) - n <ri = n (-) (-)) 

are decoherence factors with a factorized structure and for 

TO 7^ n, ruj > Uj, each factor 

,«,-(«) _ A (m, -n,) -A^/2r(m,-«,) / a 2 \ (j-, 

is expressed in terms of the associated Laguerre 
polynomials L™ {x) with the variable Aq, = 
-9j/ (2t^j)(?i-m). 

In order to show the generalized quantum thermalization, 
we need to study how the diagonal elements of approach 
the Gibbs distribution P„ cc = cxp(— /3e„), while the 
off-diagonal elements F„m vanish as a quantum decoherence 
effect. To this end, we use the basic assumptions for statistical 
mechanics that d ^ oj, but 6 ^ ujj and N ^ oo. 

First, we consider the diagonal elements. Let us establish 
the formula 

Pn [E, k) = — . . . ■ (8) 

in some sense under the reasonable assumption that 
\C{n,{nj})^ are random variables with an identical 
distribution. In fact, if this condition is satisfied, 
then these random variables have the same mathemat- 
ical expectation {Qn+i {E,6^ n)) ^ due to the restric- 
tion 1 = X)' 1^ ("-J {^j})l ■ '^hen follows from 
the law of large numbers that Y^' \C {n,{nj})\^ ap- 
proaches VIn {E ~ em5, k) {E,5^ n))^^ in proba- 
bility. Namely, the above formula is valid with high 
probability. Here we remark that with a straightfor- 
ward calculation (e.g., in Ref. this formula can 
be obtained from the microcanonical state Pmc = 
^'[1/51 AT+i [E, (5, k)] \n, {uj}) {n, {nj}\ of the "universe". 

We are now in a position to consider how P„ (Ej^) leads 
to the canonical distribution. According to Refs. |lll|2|], when 
there is no interaction, P„ {E, k = 0) does give rise to the 
canonical distribution for almost all the pure states j?/'} of the 
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constrained "universe". But when there exists an interaction 
between S and _B , is it still the case? The answer seems to 
be positive when the interaction is weak. We will attack this 
problem by theoretical analysis and numerical simulation. 

If we can show that P„ {E, k) possesses a "conformal in- 
variance" with respect to the geometrical deformation of the 
energy shell caused by the interaction, i.e., P„ (E, k) ^ 
Pn {E, K — 0), then the problem is solved. We try to jus- 
tify this "conformal invariance" for weak interaction (k ^ 0). 
It follows from the direct sum decomposition V (E, 6) = 
of the Hilbert space V {E, S) that the di- 
mension of V{E,6) can be written as Hn+i {E,6,k) = 
J2n=i ^Nin) ■ Here, ili\f{n) stands for the number of states 
in the area defined by Eq. (|3]l, and M is an upper bound of the 
summation range, which is determined by the positiveness of 
both the bath energy and the system energy. By straightfor- 
ward calculation we obtain 



{E 



(9) 



for very small S. Thus the diagonal elements of pg take the 
form 



Pn = 



[E - e„ (k)] 



N-l 



EnUE-enin)] 



N-l 



(10) 



with a finite AI and a large N. Since the eigen-energy e„ of 
the system is much smaller than the total energy of the en- 
ergy shell, we have [E - e„ (k)]^"^ « exp[(iV - 1) ln(£: + 
kA^)] exp(— /3„e„) where /9„ is the quasi-temperature defined 
as 



/3„ 



dS (E ) , 
dE 



\E^E+K.\l 



N-l 

E + k\1 



(11) 



Usually the energy correction kA„ = A„I]j \gj\ / (4wj) 
is much smaller than the shell energy E, thus the quasi- 
temperature becomes independent of 71 : /3„ w N/E. There- 
fore, Pn is "conformably invariant" with respect to k, and 
for this reason we can expect the Gibbs distribution P„ = 
exp(-/3„e„)/E^!ii exp(-/3„e„)]. 

It should be pointed out that the above argument is only 
heuristic since neither the sum ^ fli\}{n) nor the distribution 
Pn is calculated analytically. Thus we will resort to numer- 
ical simulation to support our expectation. Let us assume 
the system is a harmonic oscillator In this case, A„ — n 
and the quasi-temperature of the system is exactly /3„ — 
{N — 1) / (P — KTi^) . Due to the weakness of interaction, /3„ 
is approximately equal to N/E when n is small. We choose 
E = 0.5 and uj = 10^^ and display in Fig|2]the relation- 
ship between the distribution P„ and n for different k when 
N — 50. It is clearly shown that P„ indeed decays exponen- 
tially as Gibbs distribution as n increases. Tab. Ugives the in- 
verse temperature of the system for different coupling param- 
eters K. Here, the theoretical inverse temperature is (3 = 98. 
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FIG. 2: (Color Online) Plot of InPn as a function of the number 
n of the system, when N = 50 with different k: 5 x 10~'^(n ), 
SxlQ-^fO) f^fd 5xlO~*('A). The fitting of the Gibbs distribution 
is shown by the corresponding lines. 

These numerical results demonstrate that the state obtained 
by tracing over the bath is just the Gibbs canonical state under 
some conditions. 

TABLE I: Inverse Temperature vs Interaction 





5 X lO"** 


5 X 10"^ 


5 X 10"'' 


P 


98.94 


98.85 


98.69 



Now we consider the bath-induced decoherence effect, 
which causes the off-diagonal elements of pg to approach 
zero. When the system-bath couplings can be ignored in com- 
parison with the energy level spacing and the total energy 
of the bath, any two subspaces defined by the constraint Q 
can not overlap each other and then the decoherence factor 
vanishes. In present situation, the system-bath coupling will 
weaken this decoherence for thermalization since the positive 



terms A,, 
ings 



\n-i reduce the effective system's energy spac- 



■1+1 (k) - e„(K) = €n+l 



'^(A^+l - ^n-l) (12) 



to make them comparable to the thickness & of the energy 
shell. Then two subspaces {n, k) and {m, k) of B 
can overlap each other and the off-diagonal elements Fnm will 
not vanish. However, notice that the norm of each component 
d"^ (") jjj j-jje decoherence factor is less than unity. Thus the 

mj{m) J 

decoherence factor D may still vanish in the thermodynamic 
limit N CX3. So we can say that the factorized structure 
of D enhances the decoherence [Till , and to some extent 
compensates the negative effect of interaction in thermaliza- 
tion. 

Next, we wish to point out that in the mesoscopic case, that 
is to say, N is not large enough, novel effects may arise. In this 
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case, if the system-bath coupling is strong, there will exist fi- 
nite off-diagonal elements Fnm in the reduced density matrix. 
This means quantum coherence is introduced into the usual 
thermally equilibrium state. Such a state is called a quasi- 
thermal state. 

For a two level system with single energy spacing A, the 
quasi-thermal state can be described by the reduced density 
matrix 



Ps 



P+ F 
F* p_ 



(13) 



The diagonal elements p-t = l/(l+exp(±/3A)) approach the 
standard Gibbs distributions while the off-diagonal elements 
are non- vanishing: F — F12. We can diagonalize the above 
reduced density matrix to obtain the two effective probabil- 
ities P± (F) exactly. For small F , P± (F) can be approxi- 
mated as 

P±{F)-p±Tcoth{^)\Ff (14) 

Then the von Neumann entropy Svn — 

— J2a=± (F) In Pa [F) is approximated as 

SvN^S{E)-f3A\F\\otH^) (15) 

It is observed that due to the system-bath interaction the von 
Neumann entropy explicitly deviates from the thermodynamic 
entropy 



SiE) 



f3A 



+ 1 



ln(e' 



-0A 



1), 



which is by definition the entropy of the Gibbs equilibrium 
state Pq = diag{p+,p }. Therefore, generally von Neu- 
mann entropy does not relate to the meaningful usual notion 
of temperature. In fact, there exists no good physical notion 
of temperature for a general non-thermal-equilibrium. But the 
erasure of quantum information by thermalization indicates 
where temperature enters in this matter. 

However, only for two-level system or system with ho- 
mogeneous energy level spacing can we define an effec- 
tive temperature T^g — l/jS^g lfl4ll by the ratio r{t) = 
P+{F)/P_{F) = exp(-/3^fjA) and the level spacing A. For 
such a system interacting with a finite heat bath, even in a non 
equilibrium state, we can imagine that it is in a virtual equi- 
librium state with the effective inverse temperature 

/3,ff /?+ i^cosh2(^)coth(^) (16) 

This effective temperature in the quasi-thermal state is higher 
than the usual equilibrium temperature. Such kind of quasi- 
thermal state with a bit of quantum coherence can demonstrate 



various exotic natures in thermodynamical processes. Scully 
et al 1 15] have proposed a quantum Carnot engine in which 
the bath atoms are given some quantum coherence, which can 
increase the effective temperature of the radiation field. In 
this case, though the second law of thermodynamics is not 
violated, the quantum Carnot engine possesses some features 
that are not possible in a classical case. 



In summary, the quantum kinetic thermalization of the sys- 
tem is explored, based on an exactly solved and general model 
with weak interaction between S and B and the effect of inter- 
action is demonstrated as the deformation of the energy shell; 
Based on the model, we realized the thermalization by tracing 
over the variable of the bath B. Here decoherence is used to 
account for the disappearance of the off-diagonal elements of 
the reduced density matrix of the system in contact with a bath 
of infinitely large particle number in the weak coupling limit. 
Moreover, the novel thermodynamic effects are emphasized 
that can result from the non-vanishing off-diagonal elements 
of the reduced density matrix when the bath is mesoscopic. 
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